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We present a measurement of the pseudorapidity density of primary charged particles near mid- 
rapidity in Au+Au collisions at ^/sNN = 130 GeV as a function of the number of participating 
micleons. These results are compared to models in an attempt to discriminate between competing 
scenarios of particle production in heavy ion collisions. 
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Collisions of gold nuclei at the Relativistic Heavy Ion 
Collider (RHIC) provide a unique opportunity to study 
particle production in nuclear collisions at the highest 
available energies. In a previous publication EJ, the 
PHOBOS collaboration presented results on the energy 
dependence of the pseudorapidity density of charged par- 
ticles, dNch/dr), produced near midrapidity for central 
Au+Au collisions. It showed that this rises much faster 
with energy than in pp collisions at similar energies . 
This has been explained by the increasing role of hard 
and semi-hard processes, which are understood using per- 
turbative QCD. 

A way to control the ratio of hard to soft production at 
a fixed beam energy is to vary the impact parameter, or 
centrality, of the nuclear collisions. Soft processes, which 
produce the bulk of charged particles in pp and pA col- 
lisions, scale with the number of participating nucleons 
{Npart) in the collision (|^]. Hard processes occur in the 
interactions between individual partons in the colliding 
nucleons and are expected to scale with the number of 
nucleon-nucleon collisions, N co u. This leads to an ex- 
pected scaling of dN c h/di]\ v=0 as A x N part + B x N co u. 

Data from the WA98 experiment at CERN ||] already 
indicate possible deviations from simple Npart scaling 
even at SPS energies (^/s^at = 17.2 GeV). A stronger- 
than-linear dependence on N part is observed, well de- 
scribed by a power-law, dN c /j / 'drj\ v= o cx N part , with 



a= 1.07+ 0.04. 

Theoretical models of particle production in RHIC col- 
lisions generally fall into two classes. The first incor- 
porates the expected scaling mentioned above, using a 
Glauber model calculation @ to determine the relation- 
ship between N part and N co u as a function of impact 
parameter. The HIJING model |7j as well as calcula- 
tions by Kharzeev and Nardi (KN) [|| follow this ap- 
proach. HIJING also incorporates jet quenching and nu- 
clear shadowing which modifies the scaling, leading to a 
linear rise of the normalized multiplcity dN c h/ drj/N par t 
versus N part . KN do not include these additional effects, 
the only input parameters being the fraction of particle 
production from hard processes and the earlier PHOBOS 
result. This leads to a dependence of dN c h/drj on N part 
similar to that measured by WA98. 

The second class of calculations, based on parton sat- 
uration, predict a different dependence on the nuclear 
geometry. For example, the EKRT model ||, which in- 
corporates a geometry-dependent saturation scale, pre- 
dicts a near-constant dependence of dN c h/dri/N par t as a 
function of N part . In Ref. ||, KN also perform a calcula- 
tion based on parton saturation, including the DGLAP 
evolution of the gluon structure function. They find that 
dNch/dr]/ N par t scales as ln(Q 2 /A 2 ), where Q 2 S is the sat- 
uration momentum scale which depends on the impact 
parameter. Perhaps fortuitously, this calculation is in 
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near-perfect agreement with the other KN calculation, 
suggesting that these two physics scenarios may not be 
distinguishable, even in principle. 

In this Letter, we present the results of a measurement 
of the charged particle multiplicity per participating nu- 
cleon pair near mid-rapidity, dN ch /dr)\ \ v \<i/( ^(N part )), 
as a function of N part . For this measurement, we used a 
subset of the full PHOBOS detector, which was partially 
described in Ref. jlj. 

To measure the charged particle multiplicity, we used 
two of the three silicon detector systems implemented 
in PHOBOS Jl0|], each of which has different properties 
and thus different systematic effects on the data. The 
PHOBOS spectrometer (SPEC) used for the 2000 data 
consists of two arms, one with 16 (SPEGN) and an- 
other with 6 layers (SPECP). The first six layers of 
each sub-detector subtend — 1 < rj < 2 and A0 < 7° 
around = (SPECP) and <j> = 180° (SPECN). The 
innermost of these layers has 1 mm 2 pads while the pads 
get narrower and taller as the distance from the event 
vertex increases. The PHOBOS vertex detector (VTX) 
consists of two sets (VTXT /VTXB) of two layers which 
are located above and below the beam (z) axis. Primar- 
ily designed to measure the vertex z-position (z vtx ), the 
pads have very fine segmentation along z, but are larger 
in the x (horizontal) direction. For events with z vtx = 0, 
the detector covers A</> » ±22° around 4> = 90° and 270° 
and Arj — ±.97 around r\ = 0. 

The centrality of the collisions, from which we derive 
Npart , is primarily determined using the energy deposited 
by charged particles in two sets of paddle counters located 
at ±3.21 meters from the nominal interaction point along 
the beam axis, which subtend 3 < \r)\ < 4.5. HIJING 
simulations JTlJ suggest that, on average, the paddle sig- 
nal is monotonically related to the number of partici- 
pants, as shown in Fig. [j](a). This has been verified 
by the PHOBOS data shown in Fig. |l|(b), which shows 
the correlation between the paddle signal and the signal 
from the zero-degree calorimeters (ZDCs) which are 
located at ±18 m and measure the forward-going neutral 
spectator matter. 

As a consequence of the monotonic relationship be- 
tween the paddle signal and N part , fractions of the cross 
section as selected by the paddles correspond on aver- 
age to the same fractions of the cross section selected 
by N part . To account for the fluctuations of secondaries 
produced in the apparatus as well as of N part itself, we ac- 
tually calculate (N part ) for fractions of the cross section 
selected using a full simulation of the paddle response 
based on HIJING and GEANT. This has been done for 
each of 10 bins in the most central 45% of the total cross 
section (shown in Table ||). We find that for the PHO- 
BOS setup, ignoring all sources of fluctuations leads to 
shifts in (N part ) of less than 2%. 

A major source of experimental systematic error in the 
determination of (N pa rt) arises from uncertainty in the 



efficiency of our event selection procedure (described in 
Ref. for low- multiplicity events. We have estimated 
this by studying the frequency distribution of the num- 
ber of hit paddle counters, which is sensitive to the most 
peripheral events, and comparing the results to Monte 
Carlo simulations. By performing the procedure with 
two different models (HIJING, RQMD Q) we estimate 
a systematic error of 3% and an efficiency of 97%. Un- 
fortunately, an error of as little as 3% leads to errors on 
(Np ar t) on the order of 5% for N par t < 100. This un- 
certainty accounts for about half of the total systematic 
error on the final result described below. 

It should be noted that the Glauber model calcula- 
tion implemented in HIJING 1.35 uses a Monte Carlo 
approach. In this, nucleons are randomly distributed ac- 
cording to a Woods-Saxon distribution, and interactions 
occur with a probability proportional to the overlap of 
the Gaussian nucleon density profiles. This is very simi- 
lar to the procedure used by the PHENIX collaboration 
in a recent publication [l4| . A different approach is taken 
by KN [|],[l5| , who use a numerical integration of the nu- 
clear overlap function that should in principle give iden- 
tical results as the Monte Carlo approach. However, it is 
well known that these calculations arc done in the optical 
limit to make the integrals tractable. While this approx- 
imation is reasonable for expressing N c h as a function of 
Npart or impact parameter, it is known to be inaccurate 
for estimating the total inelastic cross section jl6]]. Thus, 
for the same fraction of cross section, it can be expected 
to give different results for (N par t) relative to a Monte 
Carlo approach. In fact, we have found that the two ap- 
proaches do disagree, and moreover, cannot be brought 
into agreement by reasonable variation of the input pa- 
rameters (e.g. radius and ctnn)- The ratio of N par t from 
HIJING ("MC") over KN ("optical") is shown in Fig. | 
as a function of N par t from HIJING. 

Following the procedure described in we have de- 
termined the charged particle multiplicity, dN c h/dr], at 
midrapidity averaged over |r/| < 1. The technique is 
based on counting "tracklets", which are three-point 
tracks consisting of two points and the measured event 
vertex. In this analysis, we have used combinations of five 
effective "sub-detectors" : layers 1 and 2 of SPECP and 
SPECN, layers 5 and 6 of SPECP and SPECN, and 
the full vertex detector. The large number of differently 
positioned detectors allowed us to control the effects of 
backgrounds and was used to check the consistency of 
our analysis technique. 

For the spectrometer, the pseudo-rapidity r\ and az- 
imuthal angle (f> of all hits in two consecutive spectrom- 
eter layers were calculated relative to the primary event 
vertex. Tracklets were then constructed by combining 
pairs of hits in both layers for which the total angular dis- 
tance D = \J Srf + 8(f) 2 satisfies the condition D < 0.015, 
where 8rj and 5<p are the deviations in pseudorapidity and 
azimuthal angle (in radians) of the two hits, respectively. 



2 



If two tracklets share a hit, the one with the larger value 
of D is discarded. A similar tracklet finding algorithm 
was used for the vertex detector. Tracklets were chosen 
in the vertex detector as combinations of hits in the two 
detector layers with \6<j>\ < 0.3, and \Srj\ < 0.04. We 
did not use the same measure D as for the spectrometer 
since the granularity in the vertex detector is substan- 
tially coarser in the <p direction. 

To study the effect of combinatorial backgrounds, we 
analyzed the full data sample with the inner layers of 
each set of detectors rotated about the beam axis by 
180 degrees. While preserving the gross features of the 
events (e.g. flow), the tracklets extracted from this data 
set arise exclusively from random coincidences of two hits 
that satisfy our quality cuts. Outside of the cut region, 
we find that the distribution of track residuals (D for the 
spectrometer, Sr) for the vertex) for the mixed-hit track- 
lets closely matches those obtained with the true detector 
geometry. By normalizing the two distributions outside 
the cuts, we thus obtain an estimate of the combinatorial 
background in the region of accepted tracklets. 

The high segmentation of the spectrometer in both 
pseudorapidity and azimuthal angle leads to a small num- 
ber of background tracklets, which is substantially larger 
for the vertex detector because of its larger pads. In 
both cases, the background level was found to scale with 
the number of occupied pads. In the spectrometer, the 
background varied from 1% to 15%, depending on oc- 
cupancy. The final number of tracklets is corrected for 
this combinatorial background with the measured dis- 
tribution, smoothed using a 2nd-order polynomial. The 
coarser segmentation of the vertex detector leads to a 
larger contribution from combinatorial tracklets. How- 
ever, since the scaling of the background with the num- 
ber of occupied pads is similar in data and simulation, we 
use a global correction factor to take this into account. 
Making an explicit correction similar to the spectrometer 
has a negligible effect (less than 1%) on the final answer. 

The proportionality factor a between the number of 
tracklets and the multiplicity of primary charged parti- 
cles for 1 77 1 < 1 was calculated using the results of simu- 
lations. Particles generated by HIJING were propagated 
through GEANT 3.21. The resulting simulated signals 
were smeared to account for detector resolution, and sub- 
jected to the same analysis chain as the real data. 

The precise manner in which these proportionality fac- 
tors were calculated and used is somewhat different in 
the spectrometer and vertex detectors. In the spectrom- 
eter, a was computed as a function of z v t x and central- 
ity. Since the spectrometer acceptance is forward of mid- 
rapidity, we only include tracklets within a fiducial cut of 
< 77 < 1. Using these proportionality factors, the cor- 
rected tracklet multiplicity was calculated using events 
in —4 < z vtx < 12 cm. For each centrality and ver- 
tex bin, the background fraction, which depends on the 
number of occupied pads, was also averaged over the se- 



lected events and then applied to the average number of 
tracklets. In the vertex detector, a was computed for 
tracklets in \r]\ < 1 as a function of z vtx (for \z v t x \ < 12 
cm) and N outeri where N outer is the number of hits in 
the outer vertex layer. This corrects for both the recon- 
struction efficiency and the combinatorial background. It 
is more difficult to distinguish these these two effects in 
the vertex detector, which lacks the pointing accuracy 
of the spectrometer. The final estimate of dNch/dr/l^^i 
was then determined in both cases by averaging the cor- 
rected number of measured tracklets over z v t x in each 
centrality bin. 

The systematic error for the spectrometer is dominated 
by the accuracy of the vertex determination and the ef- 
ficiency of the tracklet reconstruction procedure and is 
estimated to be 3%. The uncertainty on the combina- 
torial background subtraction has been estimated to be 
1%. Finally, the uncertainty in the effect of non-vertex 
backgrounds (which includes weak decays) has been es- 
timated to be less than 1%. Thus, we estimate an error 
of 4.5%, independent of centrality. 

The vertex detector has larger systematic errors, due 
to its coarser segmentation in the x-direction, which 
makes it less robust against contamination from non- 
vertex backgrounds, such as delta electrons and weak 
decays. By considering how the final result varies with 
changes in the quality cuts and event generator used, we 
estimate a final systematic error of 7.5%. 

The final value of dN c h / drj^^i is based on an aver- 
age of two separate measurements, one combining the 
four spectrometer measurements and one with the full 
vertex detector. Since the systematics are very different 
for the two different analysis techniques, we combine the 
two results weighted by the inverse of their total system- 
atic error squared to obtain the final results, which are 
shown in Table |. 

The scaled pseudorapidity density dN ~ c h / 'cM|r ( |<i / 
{.\{N P art)) as a function of N part is shown in Fig. [|, 
with (N part ) derived using HIJING. The two different 
sources of systematic error, one from the tracklet mea- 
surement and the other from the estimation of N par t, are 
combined in quadrature and shown as a band around the 
data points. The error on the participant estimation is 
based on the assumption that we may have mis-estimated 
the total cross section by 3%, our systematic error on this 
quantity. For comparison, we show an extrapolation of 
pp measurements to y^ijvjv = 130 GeV using a procedure 
described in |2| (solid circle), as well as the PHOBOS 
measurement of the pseudorapidity density at midrapid- 
ity |]lj (solid square). Our results are in good agreement 
with a recent PHENIX publication @] . 

Three model comparisons are also shown. The HI- 
JING model (solid curve) interpolates almost linearly be- 
tween the pp point and the previous PHOBOS point [0. 
Above this, we show both KN results as a single dotted 
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curve. Note that the absolute scale in this model was 
normalized to the PHOBOS value for the 6% most cen- 
tral events Finally, the saturation results of EKRT 
(dashed curve) are nearly constant as a function of N part . 
The data appear to disfavor the HIJING and EKRT re- 
sults. However, they broadly agree with the KN results, 
which are consistent with the simplest scaling expected 
by a Glauber model including a contribution proportional 
to N coU . 
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FIG. 1. a.) Simulated paddle signal as a function of the 
number of participants, b.) ZDC signal vs. paddle signal for 
PHOBOS data at yi^T = 130 GeV. 



eg . . 

St 1.1 

Q- (0 
O Q. 

? 1.05 

V 1 





30-40% 






20-30%^^- 






10-20°/^ 


-^6-10% 






^^0-6% 





100 200 300 400 

n mc 

■'part 

FIG. 2. Ratio of N part calculated by KN (optical-limit ap- 
proach) over Np ar t calculated by HIJING (MC approach) vs. 
HIJING. Each comparison is done for the same fraction of 
total cross section, as indicated next to the points. 



CM 



T PHOBOS (Glauber from HIJING) 




■ PHOBOS [1] 

• PP [2] 

EKRT [9] 

Kharzeev/Nardi [8] 
— HIJING [7] 



100 



200 



300 



N 



400 

part 



FIG. 3. The measured scaled pseudorapidity density 
dN ch /drj\ M<1 /(^(N part )) is shown as a function of N part 
(solid triangles), with N par t extracted using HIJING. The 
error band combines the error on dNch/dr)]^^ and N par t. 
The solid circle is pp data from Ref. [2]. The solid square is 
from Ref. [1]. Theoretical calculations are shown from HI- 
JING [7] (solid line), KN [8] (dotted curve) and EKRT [9] 
(dashed curve). 
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TABLE I. For each measured centrality bin, based on per- 
centile of the total cross section, we show dNch/dr/l^^x, 
the number of participants, and the final result for 
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